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ABSTRACT 

Planet formation is one explanation for the partial clearing of dust observed in the disks of some T Tauri stars. Indeed studies 
using state-of-the-art high angular resolution techniques have very recently begun to observe planetary companions in these so-called 
transitional disks. The goal of this work is to use spectra of the transitional disk object LkCa 15 obtained with X-Shooter on the Very 
Large Telescope to investigate the possibility of using spectro-astrometry to detect planetary companions to T Tauri stars. It is argued 
that an accreting planet should contribute to the total emission of accretion tracers such as Ha and therefore planetary companions 
could be detected with spectro-astrometry in the same way as it has been used to detect stellar companions to young stars. A probable 
planetary-mass companion was recently detected in the disk of LkCa 15. Therefore, it is an ideal target for this pilot study. We studied 
several key accretion lines in the wavelength range 300 nm to 2.2 pm with spectro-astrometry. While no spectro-astrometric signal 
is measured for any emission lines the accuracy achieved in the technique is used to place an upper limit on the contribution of the 
planet to the flux of the Ha-, Pay, and Pa/3 lines. The derived upper limits on the flux allows an upper limit of the mass accretion 
rate, log(M acc ) = -8.9 to -9.3 for the mass of the companion between 6 M Jup and 15 M j up , respectively, to be estimated (with some 
assumptions). 
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1. Introduction 


Transitional disks (TDs) around young stellar objects (YSOs) 
can be defined as accretion disks with an inner region that is sig¬ 
nificantly lacking in dust. Grain growth and clearing by a planet 
are two possible explanations for the removal of dusty mate¬ 
rial from an accretion disk (Alexander et al.||2013| >. Therefore, 
YSOs with TDs offer an interesting opportunity to study planet 
formation at the earliest stages (Huelamo 20131. TDs are char¬ 
acterised by a lack of mid-infrared (IR) emission and a rise into 
the far-IR (Merfn et al. 2010| . The clearing of material can be 
discussed in terms of a hole or a gap. A TD is described as hav¬ 
ing a gap when both an inner and an outer disk is present and 
a hole when there is no inner disk. Possible planetary compan¬ 
ions have very recently been identified in the TDs of two T Tauri 
stars (TTSs), T Cha and LkCa 15 ( [Huelamo et al.[[201 1 ||Kraus &| 
|Ireland|r2012| i. Both T Cha and LkCa 15 were selected for these 
studies as their spectral energy distributions (SEDs) showed no 
evidence of dust in their inner disks. The dust could have been 
cleared by a planetary mass object (|Brown et al. 2007; Espaillat 


et al. 2007 1 and these detections were done using aperture mask¬ 
ing techniques ( |Tuthill et al. |2010| i. 

Much can be learned about TDs and planet formation 
through the study of the accretion properties of YSOs with TDs. 


* Based on Observations collected with X-shooter at the Very Large 
Telescope on Cerro Paranal (Chile), operated by the European Southern 
Observatory (ESO). Program ID: 288.C-5013(A). 


For example, measurements of the mass accretion rate ( M acc ) of 
a sample of TDs by Manara et al. (2014) have shown M acc to be 
comparable to M acc in Classical TTSs (CTTSs). CTTS are Class 
II low-mass YSO with full accretion disks and strong accretion 
rates ( M acc > K) 9 M Q / yr; Muzerolle et al. 2003). The TD phase 
of a TTS is postulated to be intermediate between the CTT phase 
and the Class III phase where the disk is mostly dissipated and 
accretion has stopped. The results of Manara et al. ( 2014[ ) reveal 
information about the gas content of TDs in the regions close 
to the star and imply that these regions are gas-rich. As planet 
formation models also predict accretion onto recently formed 
planets investigations of M acc in accreting objects at the plan- 
etary/brown dwarf boundary are also important for the theory 
of planet formation (|Zhou et ah 2014). For example, Lovelace 


|et al.| ( 20TT] l explain that magnetospheric processes may govern 
accretion onto young gas giants in the isolation phase of their 
development when they have already cleared gaps in their sur¬ 
rounding disks. In this evolutionary stage they are expected to 
accrete slowly through their magnetic field lines like low-mass 
stars and it should be possible to detect strong emission in accre¬ 
tion indicators like the Ho-. [Owen | ( |2Q 14) > postulate that the high 
mass accretion rates measured for YSOs with TDs offer a clue to 
the origin of the dust poor inner regions of transitions disks. It is 
argued that the observed accretion rates in TDs imply a high ac¬ 
cretion luminosity originating from the forming planet of > 10 3 
L g . The models of [Owen (2014]) require that the planets accrete 
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Fig. 1. Full X-Shooter spectrum of LkCa 15 from the epoch 1 observations. The UVB arm is plotted in green, the VIS arm in black 
and the NIR arm in blue. Major features are marked. The spectra as shown here have not been extinction corrected. 


at least half of the material flowing into the gap in-order to trap 
sufficient sub-micron dust particles. 

Here we report on X-Shooter observations of LkCa 15. Three 
epochs of X-Shooter data were obtained between December 
2011 and March 2012 (Table 1). The primary objective of this 
project was to investigate the possibility of detecting the plan¬ 
etary companion discovered by Kraus & Ireland| ( 2012j > using 
spectro-astrometry (SA) and thus to test if SA can be used to 
detect recently formed planets (see Section 2.2). The motiva¬ 
tion to use this technique is that planets are formed very close to 
the central star and it is difficult to detect and spatially resolve 
them at optical wavelengths. This paper is divided as follows. 
Firstly, the target, the technique of SA, and the X-Shooter ob¬ 
servations are described in Section 2. In Section 3 the results 
of the spectro-astrometric analysis are outlined. In Section 4 the 
accretion properties of the star and companion are investigated 
and the results of the SA are used to estimate upper limits to the 
flux of important accretion tracers, e.g. Her, and thus upper lim¬ 
its to M acc . The conclusions are given in Section 5 and further 
information on the SA and the artefacts which we found to be 
affecting the X-Shooter data are given in the Appendix. 


2. Target, spectro-astrometry, and observations 

2.1. Target 

LkCa 15 (04 A 39 m 17L8, +22°2l'03."2) located in the Taurus- 


Auriga star forming region (d = 145 ± 15 

3C, Loinard et al. 2005) 

is a K2 TTS (Kenyon & Hartmann 

1995 

with an estimated age 

of ~ 2 Myr (Kraus & Ireland, 2012 

). Stud 

ies based on the analy- 


sis of the IR spectral energy distribution (SED) of LkCa 15 and 
sub-millimeter observations, showed it to have a TD with a mass 
of ~ 55 M j up an d a ga p of ~ 50 AU ([Andrews & Williams 2005 


Espa illat et al.| |2007[ | Andrews et al.||2011). LkCa 15 al so dis- 


plays a significant near-infrared excess (Espaillat et al. 2007 
and millimeter emission from small orbital radii ( |Andrews et al. 
[201 li ) which points towards the presence of an inner AU-sized 
disk. The mass of LkCa 15 has been estimated from the rota¬ 
tion curve of the disk at 0.97 ± 0.03 M 0 ( [Simon et al.| |2000) . 
|Manara et al.| ( |20T4j l estimate M* at 1.24 ± 0.33 M Q using the 
evolutionary models of |Baraffe et al. ( 1998| >. [Kraus & Ireland 
( 2012| > reported the detection of a possible planetary mass com¬ 
panion (LkCa 15b) at a distance of ~ 16 AU to 21 AU from 


LkCa 15, placing it at the centre of the previously detected gap. 
They estimate the mass of LkCa 15b at ~ 6 M j up to 15 M Jup . 
Specifically |Kraus & Ireland] ( |20 12} observe a blue (K or 2.1 
pm) point source surrounded by co-orbital red (L or 3.7 pm) 
emission that is resolved into two sources. They suggest that the 
most likely geometry is of a newly formed gas giant planet that 
is surrounded by dusty material, and which has been caught at 
its epoch of formation. They refer to the blue point source as 
the central source (CEN) and the red sources as the north-east 
(NE) and south-west (SW) sources. [Ireland & Kraus| ( |20T4| re¬ 
port a clockwise orbital motion of 6.0 ±1.5 degrees per year for 
the CEN source. Finally M„ rr has been estimated for LkCa 15 at 
3.6 x 1(L 9 Mq/ yr~ ( jlngleby et ak||2013| > and 4 x 10 9 M 0 / yr 
( |Manara et al.|[2UT4} . 


2.2. Spectro-astrometry 


SA is a technique by which spatial information with a high pre¬ 
cision can be recovered from a simple seeing limited spectrum 
(Whelan & Garcia 2008[ Bailey |1998| . Thus is it a very pow¬ 
erful method by which the limits placed on the spatial resolution 
of a spectroscopic observation by the seeing can be overcome. 
The technique is based on measuring the centroid of the spatial 
profile of an emission line region as a function of wavelength 
with respect to the continuum centroid. Typically this has been 
done using Gaussian fitting to produce what is referred to as a 
position spectrum. For Gaussian fitting the precision to which 
the centroid can be measured is given by 


FWHM 
2.3548( VAy 


(1) 


It is the fact that the precision of this technique is mainly re¬ 
liant on the number of detected photons (N p ) or the signal to 
noise (S/N) ratio of the observation that makes it so useful. This 
statement assumes that observations of the source under inves¬ 
tigation, where a large number of source photons are detected 
and where the detector and background noise is not significant, 
are easily achieved. Other considerations when doing a spectro- 
astrometric analysis are any curvature in the spectrum due to 
instrumental effects and the contribution of the continuum emis¬ 
sion to the results. It is straight-forward to fit the position spec¬ 
trum to remove the instrumental curvature before considering 
the final results. The continuum emission will contaminate the 
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Fig. 2. Key accretion tracers detected in the LkCa 15 X-Shooter spectrum. 


position spectrum in that any displacement in the emission line 
region will be reduced by a factor equal to (I/„„ + I cont ) / I/,,,,.. 
Studies have dealt with the issue of continuum contamination by 
subtracting the continuum emission (Whelan et al. 2004) or by 
multiplying the position spectrum of the line (after the removal 
of instrumental curvature) by this factor (Takami et al. 200l)i. 


Some of the first uses of SA were for the detection of binaries 
that were not well separated and thus were difficult to directly 
detect ( jBailey 1998| |Takami et al. 2003 i and for the detection 
of small scale jets in emission lines that were not traditionally 
believed to have a strong jet component (Whelan et al., [2004] ). 
It has since been used in further interesting ways including to 


study Keplearian rotation in circumstellar disks (Pontoppidan et 
|ah |2008|), black holes (Gnerucci et al.| 201 3[), planetary nebulae 


(B 


anco Cardenas et al. 


2014]) and jets from brown dwarfs (BDs) 


(Whelan etal. 2005 2009b). While SA only recovers spatial in¬ 
formation along the position angle (PA) used for the slit, it is 
possible to do 2D spectro-astrometry and thus assign a proper 
on-sky direction to measured offsets. Firstly, this can be done 
by taking observations using orthogonal slit PAs and combining 
the offsets measured to estimate a direction and PA for the spa¬ 
tial offsets. This method has been adopted to assign directions 
and PAs to BD outflows ( | Whelan et al.||2014||2012) . Secondly, 


2D SA can be acheived in a more direct way by analysing inte¬ 
gral field spectra (IFS). The advantage of integral field spectra is 
that it is often procured using Adaptive Optics (AO) techniques 
and so due to the much reduced FWHM of the spatial profile an 
even better spectro-astrometric precision can be achieved. For 
example, Davies et al. ( 2010) combined SA with IFS to detect 
the outflow from the massive YSO W33A with a precision of 
0.1 mas. 


Finally, the topic of spectro-astrometric artefacts should also 
be discussed. It has been shown that false spectro-astrometric 
offsets can be introduced into an observation due to a distorted 
point spread function (PSF) caused by poor tracking or uneven 
slit illumination (Brann igan et al.| |2006) . These artefacts can 
mimic the signature from a disk or jet. To check for spectro- 
astrometric artefacts one can obtain parallel spectra for exam¬ 
ple, at PAs of 0° and 180°. Any real signatures will be present 
at both PAs and their direction will be inverted between spectra 
(Brannigan et al. 2006 |Podio et ah |2008) . A further check is to 
analyse lines such as telluric absorption lines where one would 
not expect to see any signature ( |Davies et al. 2010| Whelan et 
ah] 2005). It is possi ble to model and thus remove these artefacts 
dCade Adams et al.l|2015[|Podio et al.||2008) . 
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2.3. X-Shooter observations 


The primary datasets analysed here are X-Shooter observations. 
X-Shooter, a second generation instrument on the European 
Southern Observatory’s Very Large Telescope (ESO VLT), cov¬ 
ers a large wavelength range from the ultraviolet to the near- 
IR ( |Vemet et al. 201 1| > and is proven to be extrememly useful 
for studies of accretion in YSOs ( [Alcala et al.||2014[|Whelan et| 
al. 2014a|?| l. Using X-Shooter in long-slit mode three epochs of 
data of LkCa 15 were obtained. LkCa 15 was observed in nod¬ 
ding mode with a single node exposure time of 220 s for the 
ultraviolet (UVB) and visible (VIS) arms and 200 s for the near- 
IR (NIR) arm. This yielded a nominal exposure time of almost 
30 mins after one ABBA cycle. The slit widths of the UVB, 
VIS and NIR arms were 078, 077 and 079, respectively. This 
choice of slit widths resulted in spectral resolutions of 6200, 
11000 and 5300 for each arm, respectively. The pixel scale is 
0716 for the UVB and VIS arms and 0/21 for the NIR arm. In 
each epoch two spectra were obtained. In the first epoch (El) 
spectra were obtained parallel and anti-parallel to the PA of the 
Kraus & Ireland ( 2012] ) CEN component which is assumed to be 
the PA of the companion. We take this PA to be 333°. This value 
comes from the average of the PAs reported for the K band emis¬ 
sion in Table 2 of Kraus & I reland (2012). The orbital motion 
reported in Ireland & Kraus (2014) translates to a ~ 8° change in 
the PA of the CEN source between the |Kraus & Ireland| ( |2012| > 
observations and the data presented here. Therefore Lk Cal5 b 
would be at ~ 325° at the time of the X-Shooter observations 
and therefore would still lie within the slit in our observations. 
For El the purpose of the anti-parallel spectra was to check for 
spectro-astrometric artefacts. In epochs 2 and 3 (E2, E3) spectra 
were obtained parallel and perpendicular to the PA of the com¬ 
panion. The purpose here was to recover the PA of any spectro- 
astrometric signature. Taking one ABBA cycle as one spectrum 
a total of six spectra of LkCa 15 were obtained. All this infor¬ 
mation is summarised in Table 1. 

The X-Shooter data reduction was performed indepen¬ 
dently for each arm usin g the X-Shooter pipeline version 1.3.7 
(Modiglia ni et al.[ 2010) ^The pipeline provides 2-dimensional, 
bias-subtracted, flat-field corrected, order-merged, background- 
subtracted and wavelength-calibrated spectra. Flux calibration 
was achieved within the pipeline using a response function 
derived from the spectra of flux-calibrated standard stars 
observed on the same night and of LkCa 15 . Following the 
independent flux calibration of the X-Shooter arms, the internal 
consistency of the calibration was checked by plotting together 
the three spectra extracted from the source position and visually 
examining the superposition of overlapping spectral regions at 
the edge of each arm. Any misalignment was corrected for by 
scaling the UVB and NIR spectra to the VIS continuum level. 
To correct for slit losses the spectra were normalised to the av¬ 
erage photometry available in the literature. An average V-band 
magnitude of 12.15 mag was adopted ( jGrankin et al.||2007) . The 
normalization was done by applying a correction factor equal to 
1.8. The correction factor was calculated by comparing the spec¬ 
troscopic flux, Fspec(5483 A)=2.8e-14erg/s/cm 2 /A, with the 
photometric flux F phot (5483 A)=F 0 (5483A)-10- a41215 =5.1e- 
14 erg/s/cm 2 /A, where the zero magnitude flux 


1 An improved pipeline running in a new environment i Freudling et 
|ak||2013) has been released after the conclusion of our data reduction. 
However, we have conducted tests with the X-Shooter pipeline team at 
ESO and found that these improvements do not affect the results pre¬ 
sented in this paper. 


Fq( 5483 A)=3.67e-9 erg/s/cnr/A(Iohnson 


1964 


1965). 


We have confirmed that the correction factor is the same for 


all the wavelengths where photometry is available (B & R 
from |Grankin et al.| ( |2007) l, and J, H, & K from the Two 
Micron All-Sky Survey (2MASS). Finally, the correction for 
the contribution of telluric bands was done using the telluric 
standards observed with the same instrumental set-up, and close 
in airmass to LkCa 15 . More details on data reduction, as well 
as on flux calibration and correction for telluric bands can be 


found in 

|Alcala et 

al. (2014). 



Epoch 

Slit PA (°) Date 

Time (UT) 

Seeing (") 

1 

333 

2011-12-02 

03:52 

0.85 

1 

153 

2011-12-02 

04:22 

0.84 

2 

333 

2012-02-27 

00:24 

1.3 

2 

243 

2012-02-27 

00:55 

1.5 

3 

333 

2012-03-07 

00:11 

1.15 

3 

243 

2012-03-07 

00:43 

1.35 


Table 1 . X-Shooter Observations of LkCa 15 . Each epoch con¬ 
tains two observations. In El data was taken parallel and anti¬ 
parallel to the companion PA. In E2 and E3 the slit was parallel 
and perpendicular to the companion PA. The seeing is for the 
VIS arm and it should be noted that the seeing can vary between 
the three arms. 


3. Results 

3.1. Spectro-astrometric analysis 

In Section 2 the various steps involved in the reduction of the 
X-Shooter data by the pipeline are outlined. The final product is 
a 2D spectrum where the sky has been subtracted. The spectro- 
astrometric technique is described in Section 2.2. The goal here 
was to investigate if it is feasible to detect LkCal5 b through 
the spectro-astrometric analysis of emission lines which could 
also be tracing this companion. The detection of the planet is 
theoretically possible for the same reason as it is possible to sep¬ 
arate two close stars using SA |Takami et al.112003] ). If the planet 
emits in an accretion tracer, due to accretion onto it, one will 
see an offset in the position spectrum of the accretion line at any 
point in the line where the emission from the planet is not dom¬ 
inated by the emission from the star. The size of the offset ( S ) 
will depend on the relative strength of the line flux from the star 
(Et) and planet (F p i) and the separation between the star and its 
companion (D). This assumes that the continuum emission has 
been removed (see Section 2.2). 

It was expected that the goal of detecting LkCa 15 b would 
be challenging mainly due to the fact that any contribution to the 
accretion tracers from the planet would be small. Furthemore, 
emission from transitional disk objects in accretion tracers like 
Her is known to be variable and therefore the chances of de¬ 
tecting the companion would be greatly improved if the star is 
observed in a state where ( F,) is at a minimum. Having three 
epochs of data increased the chances of observing (, F„) at a min¬ 
imum. It was found that a Voigt fit provided a better fit to the 
wings of the spatial profile than the Gaussian fit. Hence, for the 
analysis of the X-Shooter data we use a Voigt fit at all times. 
The precision in the Gaussian fit calculated using Equation 1 is 
equivalent to the 1-cr scatter in the position spectrum of the con¬ 
tinuum emission. By comparing the 1-cr scatter in the continuum 
calculated using the Gaussian and Voigt fits it was shown that the 
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Fig. 3. Spectro-astrometric analysis of the H a line in each epoch. What is shown here are the normalised intensity profiles and the 
position spectra. The size of the bin used in the analysis is given. The black line is the data along the estimated companion PA and 
the red line the data parallel and perpendicular to this PA. The dashed horizontal lines are the 1-cr uncertainty. 


precision in the Voigt fit can be well approximated by Equation 
1. Finally, it should also be noted that the spectra were binned 
when doing SA to increase the accuracy. The bin size was three 
pixels. 

The first step in our analysis was to apply SA to the UVB, 
VIS, and NIR spectra produced by the X-Shooter pipeline, in 
each epoch and at each PA. However, it was found that the UVB 
and VIS data was affected by a spatial aliasing problem which 
we found arose due to the rebinning of the spectra from the 
physical pixel space (x,y) to the virtual pixels (wavelength, slit- 
scale). Thus the processed UVB and VIS data was not suitable 
for spectro-astromertic anaysis. This problem is important for 
any future spectro-astrometric studies with X-Shooter and there¬ 
fore it is further disucssed in the Appendix. We proceeded with 
the our analysis by fitting the UVB and VIS spectra produced by 
the pipeline before the wavelength calibration step. Thus these 
spectra were bias subtracted, flat-fielded and any field cosmetics 
were removed. For the NIR arm the wavelength and flux cali¬ 
brated spectra were used. 


Aw 

Fdistort ( A X2, Ay2)^ 

Ay 

t. 





^£Fi(Axi, Ayi) 


x=0 


Fig. 4. Geometry used for the simulations of the artefacts. The 
ratio between F distort and f'i is kept constant at 1 for all the sim¬ 
ulations while Ay and </> is varied. 


that they are the lines which are least affected by absorption. 
Furthermore, the Hu line is the line in which we are most likely 
to detect the planet (Clo se et ah] 2014!. The Bry line is included 
as it is relevant to the results of|Kraus & Ireland] (|2012j>. 


3.2. Analysis of the Ha line 


We investigated the most prominent accretion tracers found 
across the three arms i.e Hy, H/i, Ha, HeIT6678, Call triplet, Hel 
1,083/rm, Pay and Pa/I. Please refer to Alcala et al. (2014) for a 
further description of these accretion tracers. In Figure [2] these 
accretion tracers are shown. The emission lines in each epoch 
and PA are presented and it is clear that the variability between 
the two PAs is not significant but is significant between the dif¬ 
ferent epochs. A study of the variability of FkCa 15 is beyond the 
scope of this paper and a second paper is in preparation on this 
subject. Position spectra were prepared for all the lines shown 
in Figure [2] and no spectro-astrometric offset was detected for 
any of these lines. The results of the analysis of the Ha and NIR 
Pay, Pa/3, Bry lines are now discussed in more depth. We chose 
to focus on the Ha, Pay and Pa/3 lines as they are the best choice 
for estimating M acc in the star and planet. This is due to the fact 


Due to the artefacts introduced by the wavelength calibration of 
the X-Shooter data the spectra before wavelength calibration was 
used for the analysis of the Ha line. Strong signatures were de¬ 
tected in all epochs and PAs and the results are shown in Figure 
[3] In El where slits parallel and anti-parallel to the companion 
PA were used the offsets do not change with PA. If this signature 
were real it would be expected that the direction of the offset 
would change between the two slit PAs (Brannigan et ah] 2006 1 . 
Therefore, what is being seen here is an artefact. Using the ap¬ 
proach of Brannigan et al. ( 2006| > the artefacts seen in the Ha line 
are mode led as being due to a distorted PSF caused by uneven 
illumination of the slit. Given a single point source with flux Fi 
located within the slit, a distortion in the PSF can be simulated 
by placing a percentage of the flux of the point source (F distort) at 
a new position off centre of the slit. This geometry is illustrated 
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Fig. 5. The spectro-astrometric artefacts detected in the raw data are simulated using equations 2 to 5. The green line is the simulated 
position spectrum. 


in Fig. [4] The centroidal position as a funcion of wavelength for 
this PSF is then given by 


f f xF(x,y)f(A + AA(y))dydy 
f f F(x, y)f(A + AA(y))dydy 


where 


F(x,y) = Gi(x,y) + Gi(x,y) 

(3) 

and 


2.773 2773 < J r-A.r 1 ) 2 +(,-Ay 1 )- 

G i(x,y) =-r-E \e <™™> 2 

y (FWHMfn 

(4) 

7 777 ry jj-j( x —A x 2)^+(y—Ay2)^ 

G 2 (x,>) - ^ FWHM)2 n F d,stone 

(5) 


In the above equations Axi, Ayi and Ax 2 , Ay 2 are the positions 
of Fi and F distort with respect to the centre of the slit. The shape 
of the signatures produced by the distortion in the PSF depend 
on the spectrum of the star while the magnitude of the false off¬ 
sets depends on the the ratio between Fi and Fdistort, the seeing 
(FWHM) and the positions Axi, Ayi and Ax 2 , Ay 2 . For their 
simulations Brannigan et al. (2006! take a value for the PA of 
the distortion ( <p ) of 45° (See Figure 4). Their results showed 
that for constant values of F i :F distort and FWHM, the magnitude 
of the false offsets increased with increasing separation (shown 
as Ay in Fig. [4]i between Fi and F distort- 

In Fig. [5] the results of the simulations of the X-Shooter Ha 
position spectra are shown. For all epochs we set Fi equal to 
F distort and it was found that in order to reproduce the measured 
artefacts it was necessary to not only vary Ay between epochs 
but also <p. Thus the PSF distortions present in our data are par¬ 
ticularly strong and vary significantly between the three epochs. 


Interestingly our results show that the size of the false offsets 
increased even as the seeing increases (see Table 1). The slit 
width remained constant at O'.'7. The full extent of the offsets 
increases from ~ 15 mas in El to ~ 55 mas in E2 and ~ 35 mas 
in E3. Brannigan et al. ( 2006} demonstrated that the size of the 
artefacts could be minimsed by increasing the FWHM with re¬ 
spect to the slit width (while keeping the parameters describing 
the PSF constant). They suggest that using a slitwidth which is 
narrow in comparison to the seeing could offer some protection 
against these artefacts. However, our results for E2 and E3 show 
that this will not solve the problem with artefacts in cases of a 
strong distortion of the PSF. 

To check for a residual spectro-astrometric signature we sub¬ 
tracted the artefacts firstly by subtracting the El anti-parallel 
spectra and secondly by using the model results. No real spectro- 
astrometric signature was seen. Although no spectro-astrometric 
signature from the planet is detected the uncertainty on the anal¬ 
ysis can be used to recover an upper limit for the flux of the Ha 


emission from the planet. Garcia et al. ( 1999| discuss how the 
spectra of two unresolved companions can be recovered using 
the following equations 


F, = (1/(1 +r) + S(A)/D)F total (A) 


( 6 ) 


F P i = (r/( 1 +r) + S 0A)/D)F tolal (A ) (7) 

where r is the flux ratio of the two components in the contin¬ 
uum, S(T) is the centroid and D is their separation. As we detect 
no spectro-astrometric signature S(T) is replaced by the 3-cr un¬ 
certainty measured for each value of A using Equation 1. To esti¬ 
mate the upper limit for the Ha integrated flux we fit the Ha line 
in the F p i spectrum. As we can only substitute the uncertainty 
for S(/t), we cannot infer any information about the shape of the 
Ha line profile in the spectrum of the companion. Therefore, it 
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is assumed that shape of the star’s and companion’s Her line pro¬ 
file is the same. This is s ufficient for estimating a n upper limit 
for M acc onto the planet. [Kraus & Ireland (2012) estimate the 


mass of LkCa 15b to be between 6 M j up and 15 M j UD . It was 


*Jup‘ 


found that for a planet with this mass range the r/(l+r) term was 
negligible compared to the S(/l)/D term. For example, consider¬ 
ing the R and I band we find that our estimate lies between 10 4 
and 10~ 2 while S(d) / D is on the order of 10 _1 . Please refer to 
Appendix B for more details on how r was estimated. Thus the 
upper limit on the flux from the planet was dependent on the 
uncertainty in the spectro-astrometric analysis and the distance 
between the star and the planet taken at 071. The 3-cr upper limit 
of the Her flux from the companion (averaged from all 6 spectra) 
is 1.0 x 10~ B erg/s/cm 2 when a correction for A,. = 1.2 (Manara 


et al. 


the companion in Section 4. 


3.3. Analysis of NIR lines 


Zhou et; 


al. (2014) carried out a study of M acc onto 


For example, 

three accreting objects, GSC 06214-00210b (GSC), GQ Lup b 
(GQ) and DH Tau b (DH). The mass range covered by these 
objects is 10 M. Jup to 30 M Jup and they are among only a very 
small number of BDs / planets for which M acc has been esti¬ 
mated. 


Zhou et al. 


(2014) estimated M acc using UV and optical 


photometry done with the Hubble Space Telescope (HST) and 
their method which follows Gullbring et al. ( 1998) assumes that 
their objects accrete material from a disk in the same manner as 
pre-main sequence stars. They found their measured values of 
M acc to be higher than expected when compared to the correla¬ 
tion between M* and M acc , derived for YSOs and BDs, from UV 
excess measurements and other conventional methods for mea¬ 
suring M acc (Alcala et al. 


2014 


Herczeg & Hillenbrand 


2008). 


This correlation would apply to objects formed like stars (proto- 
stellar core fragmentation) and these results of Zhou et al.|( |2014| ) 
could point to an alternative formation mechanism for GSC, GQ 
and DH. The authors briefly discuss formation by core accretion 
and disk instabilities. Furthermore, |Zhou et aI7|(|2014]) conclude 


that these high values of M acc are consistent with the presence of 
massive disks around these objects from which they accrete, as 
their separations from their primaries (100-300 au) are too large 
to support accretion at such high rates directly from the primary. 

We wish to add to the work of Zhou et al. ( 2014) ) by deriving 
M acc f° r LkCal5 b. To this end we use the upper limit on the 
flux from the Ha, Pay and Pa /3 lines to estimate the accretion 
luminosity ( L acc ) in each line. The equation M acc = 1.25 ( L acc 
R„) / (G M„) is used (|Gullbring et al. 1998 Hartmann et al. 


2014) is applied. This value is used to estimate M acc for 


In Section 3.4 of their paper Kraus & Ireland ( 2012) outline re¬ 
jected alternative explanations for their observations. They dis¬ 
cuss how SA of the Bry line could be used to investigate the 
possiblility that their observations are not of continuum emis¬ 
sion from a planetary companion and circumplanetary dust but, 
are of line emission from gas in the disk gap. The Bry, Pay, and 
Pa (3 emission in each epoch is shown in Fig. [2] and no spectro- 
astrometric signature is detected for any of these lines. Also note 
how faint the Bry line is and the average S/N is found to be 
6. This supports the arguments of Kraus & Ireland ( |2012| ). For 
the spectro-astrometric analysis the pipeline processed data was 
used as the NIR arm does not suffer from the spatial aliasing 
effects described in Section 3 and in the Appendix. Also arte¬ 
facts such as those seen in the Ha line were not detected. Using 
the uncertainty in the spectro-astrometric analysis of the Pay and 
Pa p lines and the method outlined above the 3-cr upper limit for 
the Pay and Pa (3 flux from the planetary companion is estimated 
at 7.0 x 10~ 15 erg/s/cm 2 and 8.0 x 10~ 15 erg/s/cm 2 , respectively. 
These values have been corrected for extinction. This flux is also 
used in Section 4 to place an upper limit on M acc for the planet. 


4. Mass accretion in LkCa 15b 

Studies of mass accretion and M acc onto planetary mass objects 
should help to distinguish their possible formation mechanisms. 


and the results of |Alcala et a l. (2014) are used to convert 
the luminosity of each line (Ln„ e ) to L acc . R p / is taken as 0.22 R 0 
and 0.34 R 0 for masses of 6 M j up and 15 M j up respectively. R p i 
was estimated using the SETTL evolutionary models for lMyr 
(Allard et al. 2003) . Our method is analogous to studies which 
use various accretion tracers to derive M acc in low mass stars 
and BDs ( |Alcala et al.| |2014[ [Whelan et al.) |2014a) , and there¬ 
fore it is assumed here that the correlation between and M acc 
for low mass pre-main sequence stars is applicable for planetary 
mass objects. Studies of M acc using L acc estimated from L/,„ e 
have shown that the optimum way to estimate M acc is to use 
several accretion tracers fou nd in different wavelength regimes 
(Rigliaco et al. 2011 2012). In this way any spread in M acc due 
to different tracers probing different regimes of accretion or hav¬ 
ing a jet/wind component can be overcome. This is our reason 
for using more than one accretion tracer here. However, as many 
of the accretion tracers found in the spectrum of LkCa 15 show 
strong variability (refer to Fig. [2j we chose to limit our calcula¬ 
tions to the lines with the least variability namely Ha, Pay and 
Pa/i Here we are specifically referring to the presence of vari¬ 
able absorption features. For example in El the He I 1.083 n m 
line is double peaked but in E2 and E3 it is single peaked with 
a deep red-shifted absorption. The Hy and H/l lines also show 
red-shifted absorption in E2. Also Ha is the line in which we 
are most likely to see a contribution from the planetary compan¬ 
ion Close et al. ( 2014) . It is found that for the 3-cr upper limit 
on the flux, log(M flcc [Ha]) = -9.0 to -9.2, log(M acc [Pay]) = -9.1 
to -9.3 and log(;V/ flCT [Pa/l]) = -8.9 to -9.1. The range in M acc for 
the individual lines comes from the estimated range in the mass 
of the companion. These estimates are in agreement with predic¬ 
tions for mass accretion rates from circumplanetary disks, which 
place M acc for a 1-10 M j up planet in the range 10~ 9 to 10~ 8 M 0 / 
yr ( jZhU][20T5i . 

To compare the estimated accretion rate of LkCal5b with 
that of the central object, we have also derived M acc for LkCal5. 
In fact, Manara et al. (2014) used the El spectra of LkCa 15 
presented here, in their study of M acc in YSOs with TDs. They 
make their estimates of M acc by modeling the excess emission 
produced by disk accretion using a set of isothermal hydrogen 
slab emission spectra. As a check on their values of the accretion 
luminosity (L acc ) from the models, they also derive L acc using 
several different accretion tracers and the L acc , line luminosity 
(L u ne ) relationships of Alcala et al. ( |2014) . Manara et al. ( |2014) 
report M acc = 4 x 10 -y M 0 / yr for LkCa 15. Here we estimate 
Macc for LkCa 15 in the same manner as described above for 
LkCa 15b and again using the Ha, Pay and Pa (3 lines. For LkCa 
15 L ij n e is corrected for extinction using A t , = 1.2 and R* and 
M, are taken as 1.52 IQ, and 1.24 M 0 respectively (Manara et 
[ak| |2014) . Our motivation is to check for any strong variability 
in M acc for LkCa 15 between our three epochs of data. In Fig. 
|6]log(M acc ) calculated using the Ha, Pay and Pa/1 lines from all 
six X-Shooter spectra are plotted. The dashed lines are the mean 
value of log(M flcc ) in each line. Considering the three lines the 
mean value of M acc is (1.3 + 0.6) x 10 -9 M 0 / yr. This is the 
same order of magnitude but somewhat lower than the values re- 
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ported in Ingleby et al.| ( [20TT) and Manara et al.| ( |20T4] >. Overall 
the variation in M acc is not significant and is within the values 
reported by Costigan et al. (|2012|> (0.37 dex) and Venuti et al. 


( |2014| > (0.5 dex) in their studies of accretion variability in large 
samples of YSOs. 

To put all these results into context we plot in Fig. 0 M acc 
derived from X-Shooter data for a number of YSOs and BDs, as 
a function of mass. The results for the YSOs come from lAlcala l 


et al. (2014 1 and the linear correlation marked by the black line 
comes from the fit to the Alcala et al. 2014[ > results. Also shown 
are the results from |Zhou et al.| ( 2014 > and |Close et al.| ( |2014| l. 
Close et al. ( 2014[ > use the Magellan Adaptive Optics system 
(MagAO) to detect a 0.25 L 0 companion in the disk of the YSO 
HD142527. They estimate M acc from the Ho emission to be 5.9 
x 10 -10 M 0 / yr and discuss how this techniques could be used to 
detect planets in their gas accretion phase. The names of the BDs 
and planets are marked on the figure. The estimates for M acc in 
LkCa 15 and LkCa 15b are marked on Fig. [7] The upper lim¬ 
its we show here for LkCa 15b are comparable to M acc for GQ 
and do not contradict the conclusions of |Zhou et al.| ( |2014] > that 
some planetary mass companions exhibit high mass accretion 
rates when compared to BDs and low mass YSOs. 


The central assumption of our study and that of Zhou et al. 

( 2014| > is that magnetospheric accretion can also be used to de¬ 
scribe accretion from a circumplanetary disk. However, it should 
be noted that it is not yet understood if magnetospheric accretion 
is the dominant mechanism by which planets accrete material 
from a disk. In their paper on observational signatures of ac¬ 
creting circumplanetary disks |Zhu ( 2015) > discuss a number of 
caveats to this assumption. If M acc is high and the central ac¬ 
creting object has a weak magnetic field, accretion will occur 
via a boundary layer rather than along magnetic field lines. In 
this scenario strong optical and UV emission from the boundary 
layer could be expected. Furthermore, if magnetic fields are high 
enough to truncate the inner disk and for magnetospheric accre¬ 
tion to occur, it is not clear if heating mechanisms are sufficient 
to produce strong line emission such as what is seen in pre-main 
sequence stars. As more and more planetary mass objects are de¬ 
tected a clearer understanding M acc at planetary masses should 
emerge. Future studies of circumplanetary disks with the James 
Webb Space Telescope (JWST) for example, should also answer 
many questions. 



Spectrum 


Fig. 6. log(M afc ) calculated for all six spectra using the Ha, Pay 
and Pa/J lines. The dashed lines are the mean values from each 
tracer. 


5. Summary and conclusions 


Here multi-wavelength spectroscopic observations of the transi¬ 
tional object LkCa 15 obtained with X-Shooter are presented. 
The goal was to investigate the possibility of using spectro- 
astrometry to detect planetary mass objects in transitional disks. 
It is planned that the data presented here will be part of an exten¬ 
sive examination of the variability of LkCa 15 and further data 
will be obtained in winter 2014-2015. Our results and concul- 
sions can be summarised as follows. 

X-Shooter is not the best choice for a spectro-astrometric 
study due to the artefacts which can be introduced by the pro¬ 
cessing of the data and by instrumental distortions of the PSF. 
The artefacts were primarily a problem for the VIS and UVB 
arms and for the broadest brightest lines, i.e Ha. X-Shooter is a 
complex instrument with flexures both from thermal and gravity 
effects (due to its Cassegrain position), in addition to the curved 
and tilted spectral format (which necessitates elaborate data pro¬ 
cessing). Instruments which do not require such a high level of 
data processing to produce the final spectrum would be a better 
choice for spectro-astrometry. 

Although no spectro-astrometric signatures were detected 
for any of the accretion lines the uncertainty in the analysis could 
be used to put an upper limit on the flux from the planetary com¬ 
panion, in each line investigated. This was done for the Ha, Pay 
and Pa/1 lines. The derived fluxes when corrected for extinction 
were used to place an upper limit on the planetary companions 
mass accretion rate. It was found that log(M acc ) = -8.9 to -9.3 
for the mass of the companion between 6 M j up and 15 M j llp , re¬ 
spectively. This is in agreement with estimates of mass accretion 
rates from circumplanetary disks. However, a more accurate es¬ 
timate is needed before any conclusions can be drawn about the 
formation of LkCa 15b. 

Overall we conclude that spectro-astrometry may not be the 
most efficient mechanism for detecting young planets in tran¬ 
sitional disks. Even without the difficulties caused by spectro- 
astrometric artefacts, observations would have to be made when 
Ft / F p i is favorable for a detection and sub-milliarcsecond accu¬ 
racy would be needed. However, the results we present here do 
show that it should be possible to detect LkCa 15b in Ha with 
high angular resolution imaging techniques for example, with 
VLT / SPHERE or MagAO. Taking the maximum value of the 
Ha flux from the star as 1.54 x 10 12 erg/s/cm 2 (measured for 
the E3 333° spectrum), the upper limit on the Ha flux from the 
companion would give a contrast of ~ 3.4 mag between the star 
and planet. It is feasible to detect this with SPHERE, which pro¬ 
vides a contrast of ~ 6.5 mag at 80 mas in 1 hour of exposure 
time, in the Ha narrow band filter. This is still feasible even if 
we take the flux from the planet as being 50 times less than the 
upper limit reported here. These observations could then be used 
to better estimate M acc . Furthermore, MagAO has already been 
able to spatially resolve a companion in Ha within the gap of 
the transitional disk of the Herbig Ae/Be star HD 142527 (Close 
et al. 2014 1 . Here Amag = 6.33 ± 0.20 mag at the Ha line. A 
detection of LkCa 15b in Ha would allow the estimate for M acc 
given here to be greatly refined. 
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Fig. 7. Comparing the upper limits for M acc in LkCa 15b with values derived for accreting objects with masses within the range of 
low mass stars to planets. All the BDs and planets are marked with their names. For the objects from the Zhou et al. ( 2014| ) study 
M acc was estimated using UV and optical photometry. For the Close et al. (2014) object M acc was estimated from the Ha emission 
detected with MagAO. For the rest of the objects M acc was derived from X-Shooter data using the method outlined in Section 4. 
The majority of the sources are low mass stars and the values of M acc are taken from Alcala et al. (2014). The black line represents 


the linear fit to the results of 
1-cr deviation from the fit. 


Alcala et al. 


(2014) and it shows the correlation between M and M acc . The dashed lines represent the 
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The above ratio was estimated assuming only photospheric emission, that 
is, not considering accretion. Firstly we stress that the contrast has been calcu¬ 
lated in the R-band where the continuum excess due to accretion onto a planetary 
mass object is much less than in the blue (see Fig.2 byjZhou et al. ( 2014)). From 
the slab models shown in Fig 2 ofJZhou et al. 12 014) one can see that the con¬ 
tinuum excess in the R-band is in most cases just a fraction of the photospheric 
emission. The exception is GQ Lup where the excess continuum emission is 
several times the photospheric emission at a wavelength close to the Ho line. 
Taking the example of GQ Lup as the extreme case we place a factor of 10 as 
the upper limit on the ratio between the continuum excess emission and the 
photospheric emission. Increasing r by a factor of 10 increases the range we 
estimate for r to between ~10 2 and 8x10 4 . Thus r is still negligible com¬ 
pared to all values of S(/l) / D, and we find that this increase does not change 
our estimate for the upper limit for the Ho line flux from the planet. The 
continuum excess emission due to accretion onto the planet, is also negligible 
in the NIR (see again Fig.2 in |Zhou et al. 1)2014) ). Therefore, the contrast is not 
affected. 


Appendix A: Spectro-astrometric analysis of the 
pipeline processed data. 

The full UVB, VIS and NIR in each epoch and at each PA were analysed using 
spectro-astrometry. In the VIS arm and the brightest parts of the UVB, an oscil¬ 
lating pattern with an amplitude of ~ 10 mas was seen. To illustrate this better 
the position spectra for all six spectra (three epochs x 2 PAs) for the regions of 
the H/5 (panel a), Ha (panel b) and Pa/5 lines (panel c) are presented in Fig. |A.l| It 
is concluded that this signal is caused by spatial aliasing introduced by the data 
reduction. This signal is not present in the raw data (refer to Fig.[3j. Four to five 
pixels are normally necessary to adequately sample the FWHM of a spatial pro¬ 
file given by the seeing. If under sampled the effects of spatial “aliasing” appear, 
which show-up as spatial oscillations in the ID spectra and in their correspond¬ 
ing 2D frames. The problem arises here due to the rebinning of the spectra from 
the physical pixel space (x,y) to the virtual pixels (wavelength, slit-scale) and 
is enhanced for high S/N and in the wavelength range of highest instrumental 
sensitivity, which in this case is the VIS arm. The sensitivity to S/N is the reason 
why we do not see this effect in the NIR arm and only in the brightest parts of 
the UVB arm. Also this effect will be enhanced in good seeing conditions. For 
example, note that in the UVB arm the effects are worse in the El spectra which 
has considerably better seeing than the E2 or E3 data (see Table 1). 

From Fig. |A.1| it is clear that this effect will have important consequences 
for any spectro-astrometric analysis. Being intrinsic to the interpolation method 
along the slit, the problem cannot be totally removed, however its effect can be 
reduced. This was done by using a smaller bin than the standard one to interpo¬ 
late the pixels along the slit. Therefore, in order to minimize the aliasing effect 
on the spectro-astrometric analysis, we have used a value of O'/05 as input for the 
binning parameter “rectify-bin-slit” instead of the standard value of O'/16. As is 
seen in Fig. |A. 1| panel d, this approach solved the spatial aliasing problem in the 
VIS arm. However, it should be considered that the re-binning in order to solve 
the aliasing problem could alter any real spectro-astrometric signal. 


Appendix B: Continuum flux ratio of the LkCa 15 
it’s companion 

To estimate the continuum flux ratio of LkCa 15 and it’s planetary companion 
(r) we need to derive the flux of the star and planet at the given wavelength. In 
the case of the star and the Hq. line, we have considered the faintest Johnson 
R-magnitude from Grankin et al. (2007), since it represents the quiescent state 
of LkCa 15 (R= 11.97 mag, or M/?= 6.23 mag for a distance of 140 pc). In the 
case of the planet, we have estimated for a 15 and a 6Mj up object using 
the SETTL evolutionary tracks (Johnson system) for 1 Myr. We obtain values of 
13.486 and 16.47 mag, respectively, and a difference in magnitudes of 7.256 mag 
and 10.24 mag with respect to the star. We estimate r values of ~ 10 -3 and 8x 10 -5 
for 15 and 6M j up . Therefore, for a planet in this mass range, ther/(l + r ) term 
is negligible. For the Paschen lines, we have repeated the same procedure using 
the 2MASS J-mag of the star (J=9.42±0.02 mag), and the 1 Myr SETTL models 
in the 2MASS system to estimate the magnitude of the planet. We estimate a 
difference in magnitudes, A J, between the planet and the star of 5.09 mag and 
6.845 mag for a 15 and a 6 M j up object, respectively. The r values are ~ 9x 10 -3 
and ~2.10 -3 , so that the r/( 1 + r) term is also negligible. 
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Fig. A.l. Spectro-astrometric analysis of a section of the UVB, VIS and NIR spectra. Note the small-scale oscillations in the UVB 
and VIS spectra. This is caused by spatial aliasing and the effect is minimised by reducing the bin used to interpolate the pixels 
along the slit. 
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